Probing Neutral Majorana Fermion Edge Modes with Charge Transport 
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We propose two experiments to probe the Majorana fermion edge states that occur at a junction 
between a superconductor and a magnet deposited on the surface of a topological insulator. Combin- 
ing two Majorana fermions into a single Dirac fermion on a magnetic domain wall allows the neutral 
Majorana fermions to be probed with charge transport. We will discuss a novel interferometer for 
Majorana fermions, which probes their Z2 phase. This setup also allows the transmission of neutral 
Majorana fermions through a point contact to be measured. We introduce a point contact formed by 
a superconducting junction and show that its transmission can be controlled by the phase difference 
across the junction. We discuss the feasibility of these experiments using the recently discovered 
topological insulator Bi2Se3. 

PACS numbers: 71.10.Pm, 74.45. +c, 03.67.Lx, 74.90.+n 



Majorana fermions have attracted interest in con- 
densed matter physics because their exotic non-Abelian 
quantum statistics[l| form the basis for topological quan- 
tum computation [J, HI- Potential electronic systems 
hosting Majorana fermions include the v = 5/2 quantum 
Hall state [HQ, the p-wave superconductor S^RuOj 
and topological insulator/superconductor structures 
8]. In the v = 5/2 quantum Hall state, a Majorana bound 
state is associated with the charge e/4 quasiparticle, and 
gapless chiral Majorana fermions form the neutral sec- 
tor of the edge states. Thanks to the e/4 charge, the 
quasiparticle's non-Abelian statistics can be probe d by 
measuring charge transport of the edge states i| nj U | • 
Recent experiments have shown evidence for the quasi- 
particle charge e/4p3, [l3| . and there are now intense 
efforts to prove or disprove their non-Abelian nature. 

Detecting Majorana fermions in superconductors is 
more challenging because they are electrically neutral. 
In this work, we propose two experiments to probe neu- 
tral Majorana fermion edge states predicted in supercon- 
ductor/magnet/topological insulator structures [6[. Our 
basic setup, shown in Fig. 1, involves a grounded su- 
perconductor surrounded by two magnets with opposite 
out-of-plane magnetization, which are both deposited on 
the surface of a topological insulator. The magnetic do- 
main wall gives rise to chiral Dirac fermions that play 
the role of "leads" connecting the superconductor to the 
source and drain. An electron incident from the source 
splits into two Majorana fermions which take different 
paths around the edge of the superconductor and then 
recombine before going to the drain. We will show the 
source to drain conductance probes the interference of 
the Majorana fermions, forming a novel "Z2 interferom- 
eter" . In addition, we will show that the transmission of 
Majorana fermions through a "point contact" formed by 
a Josephson junction between two superconductors can 
be measured, and that the transmission can be tuned by 
controlling the phase difference across the junction. 

A topological insulator 14 , HU has gapless surface 
states that are topologically protected in the absence of 



FIG. 1: An interferometer for Majorana fermions. Magnetic 
(M) and superconducting (SC) materials are deposited on a 
topological insulator. Chiral Majorana fermion edge states 
(denoted by a single arrow) circle the outer boundary of the 
superconductor, and chiral Dirac fermion edge states (denoted 
by the double arrow) are confined to the magnetic domain 
wall connected to a source (S) and drain (D). A return path 
between the drain and source is shown with the dashed line. 
When a voltage is applied to the source electrons are split into 
two Majorana fermions, allowing their Z2 interference phase 
±1 to be probed by measuring the current in the drain. 



time reversal or gauge symmetry breaking fields. Break- 
ing time reversal symmetry either by an applied magnetic 
field or by depositing a magnetic material can open an 
energy gap leading to a novel surface quantum Hall ef- 
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Depositing a super- 
conductor on the surface leads, via the proximity effect, 
to a surface superconducting state that hosts Majorana 
fermions In view of the recent experimental discover- 
ies of topological insulator phases in Bi^Sbi-a; ll| lj| and 
Bi^Se^ [20| . and the earlier experimental evidence of good 



contact between superconducting Nb and BLSbi_-r[2ll|. 
the experimental study of these novel gapped phases is 
now possible. 

The superconducting and magnetic phases of the sur- 
face states, as well as the gapless states at interfaces 
between them, can be described with the Bogoliubov 
de Gennes (BdG) formalism. The Hamiltonian is H = 
^WS>/2, where * = (tp^, ip±, ip] , -i/jI) t and 
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Here ijj^ and ipi are electron operators of the surface 
states which are Kramers degenerate at k = 0. The 
first line in Tt describes the free surface states coupled to 
the vector potential A. a — (a x ,a y ) are Pauli matrices, 
vf is the Fermi velocity and /x is the chemical potential. 
A^|i/;j + h.c describes the superconducting proximity ef- 
fect. Spatially uniform A gives a gapped excitation spec- 
trum = y/ (±ujkj — p) 2 + A 2 . M^a z ip describes the 
Zeeman splitting due to the magnet. Spatially uniform 
M gives E£ = a/ v 2 |k| 2 + M 2 ± /x, which is gapped when 
M > /x. The BdG Hamiltonian has particle-hole sym- 
metry, expressed by = where the particle-hole 
operator is S£ = a v T v ^*. The eigenstates £±e with en- 
ergy ±E obey = and only the E > half of 
the spectrum represents independent excitations. 

An interface between two half planes (y > and y < 0) 
with different mass terms gives rise to gapless ID domain- 
wall states. First consider a superconductor-magnet in- 
terface modeled by A = A 6(y) and M = Mq6(— y). 
Solving Q, we find one chiral branch of bound states 
with a four component wavefunction £,k(x,y) localized 
near y = 0. £fe=0 has zero energy and satisfies S£o = £,o, 
which fixes its phase up to a ± sign. Using k ■ p theory 
the eigenstates for small k are £,k{x,y) — enp(ikx)£ > o(y) 
with energy E(k) = hv M k, where v M = v F (£a\T z Vy\£o) = 
v F yfl- p 2 /M 2 /{\ + n 2 /Al). These define Bogoliubov 
operators 7a, = J dxdy£k(x, y)^(x, y) which satisfy 7^ = 
7_fe. The continuum operators 7(2:) ~ J dk^e %kx are 
Majorana fields, J'(x) — j(x) obeying the low energy 
Hamiltonian H = —ihvM"fd x "f. 

To model a magnetic domain wall we take M = 
M sgn(y). We find a gapless branch of chiral edge states 
between a xy = ±e 2 /2h. When expressed in the BdG for- 
malism, two chiral branches of bound states with energy 
E(k) ~ hvjjk appear due to the double counting. For 
E(k) > 0, the two states have the form fa ® \t z = 1) 
and Qf-k ® \t z = —1). where fk(x,y) is a two compo- 
nent wavefunction in the a z sector and 0/ = <r y f* is the 
time reversal operator. These correspond to the electron 
operators cj. and C-k respectively. 

To analyze the device in Fig. 1, we employ the BCS 
mean field theory to calculate the transport current due 
to quasiparticles. This is justified because the supercon- 
ducting order parameter at the surface inherits its phase 
from the bulk 3D superconductor, which behaves classi- 
cally at low temperature. When the source is biased at 
a subgap voltage V <C Ao the quasiparticles involved are 
exclusively the gapless Majorana fermion edge states. 

An electron incident from the source can be transmit- 
ted to the drain as an electron, or converted to a hole by 
an Andreev process in which charge 2e is absorbed into 
the superconducting condensate. Before solving the gen- 
eral source to drain transmission problem we will show 
that the behavior at E — follows from a simple ar- 
gument. Scattering at the left tri-junction, where the 



incident Dirac fermion meets the superconductor, must 
transform an incident E = electron into a fermion 
if> built from the Majorana operators 71 and 72. The 
arbitrary sign of 71,2 allows us to choose ip = 71 + 172. 
Likewise, scattering at the right tri-junction transforms 
ip into a fermion in the right lead. This must be either 
c^j or cr. A superposition of the two is not allowed be- 
cause it is not a fermion operator. To determine which 
occurs, we observe that when the size of the superconduc- 
tor shrinks continuously to zero, the left and right lead 
seamlessly connect to each other. Adiabatic continuity 
thus dictates that an incident E = electron is trans- 
mitted as an electron, c\ — > c* R . When the ring encloses 
a quantized flux <f> = nh/2e, this adiabatic argument 
breaks down. Instead odd n introduces a branch cut for 
one of the Majorana modes, i.e. 71 — > —71. Thus, when 
the ring encloses an odd number of flux quanta, c\ — > cr, 
and an incident E = electron is converted to a hole. 

To obtain the scattering probabilities at finite en- 
ergy < E -C A, we use the BdG formalism to solve 
the scattering problem in the limit that the size of the 
ring L is much larger than the decay length of the 
Majorana edge states into the bulk, which is of order 
max(hvF / Ao, Hvf/Mq). First consider the scattering at 
the left tri-junction. A 2 x 2 scattering matrix S(E) re- 
lates the two incoming states in the left lead \t z = ±1), 
which we denote e and h (for electron and a hole), to 
the two outgoing Majorana edge states £1 and £2 on the 
top and bottom of the ring, (£i,6) T = S(E)(e, h) T . To 
simplify the notation, we have used the channel label to 
denote the amplitude of the scattering states in the cor- 
responding channel. Particle-hole symmetry implies that 
S(E) = S*(-E)t x . At E = 0, this property, along with 
unitary S'S= 1, allows S to be chosen as 

so that c 1 — > (71 — 172) /V2- Another solvable limit 
is when the BdG Hamiltonian has a mirror symmetry 
H{-y) = M~ 1 H(y)M with M = ia y . The electron and 
hole channels are eigenstates of A4 with eigenvalue ±i, 
whereas the two Majorana fermion edge states are in- 
terchanged. This leads to @ at any energy. To obtain 
the exact scattering matrix at E 7^ for a tri-junction 
without mirror symmetry requires solving the 2D scat- 
tering problem. Here we assume that Eq.([2]) is a good 
approximation of the scattering matrix at low energy. 

Next we study the propagation of the chiral Majorana 
fermion. When there is no magnetic flux, in the semi- 
classical limit the wavefunction at energy E <C Ao can 
be approximated by £{l,s) — £o( s ) exp(ik(E)l), where 
I parameterizes the length along the interface, s pa- 
rameterizes the distance perpendicular to the interface 
and k(E) = E/vm- In the presence of a magnetic 
flux $ = nh/2e, the superconducting phase <fi winds by 
2nn around the ring accompanied by a vector potential 
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A = V</>. It is convenient to choose a gauge in which 
the spatial variation of <f> is concentrated near the mid- 
dle of the upper semi-circle. Away from this "scattering 
region" , the wavefunction is a free Majorana edge mode 
as before. This scattering problem can be solved with 
a U(2) gauge transformation that eliminates the spatial 
variation of <f> and the nonzero A. The wavefunction of 
the exact scattering state is simply the undisturbed wave- 
function multiplied by exp[iT z </>(/)/2]. We conclude that 
in the presence of a magnetic flux $ = nh/2e the chi- 
ral Majorana edge mode 71 acquires an additional phase 
shift nir across the junction. 

The scattering amplitude of the ring is found by com- 
posing the scattering matrices: 
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The current in the drain when the source is biased at volt- 
age V and the superconductor and drain are grounded is 



/=(- 



(-1)"- / dE[f(E-eV)-f(E + eV)} cos 6(E), 

(4) 

where / is the Fermi-Dirac distribution function and 9 = 
k(h — I2) = ESL/vm is the relative phase between two 
paths of different lengths. Evaluating the integral we find 



n e irkBT sm(eV 5L / vm) 
h smh(TT ksT5L/vM) 



k B T, eV < A . 



(5) 

At fixed bias, the current "oscillates" as a function of 
the discrete magnetic flux nh/2e, reflecting the Aharonov 
Bohm phase for Majorana fermions, which takes val- 
ues ±1. Our device thus functions as a "Z 2 inter- 
ferometer" for Majorana fermions. The "visibility" of 
these oscillations is suppressed below a temperature scale 
ksTsh = Tivm/SL due to thermal averaging. In addition, 
at finite bias voltage the current oscillates as a function 
of V with a period 2irkBT$L/ 'e due to the energy depen- 
dence of the relative phase. That the oscillation persists 
to high bias voltages without any damping is due to the 
absence of dephasing in our calculation. A similar sit- 
uation occurs in the electronic Mach-Zehnder interfer- 
ometer: the decay of the magnitude of interference os- 
cillation at high bias voltage is attributed to dephasing 
processes [2^|. Sources of dephasing in our system include 
coupling of Majorana fermions with other degrees of free- 
dom, as well as interactions between Majorana fermions. 
Since Majorana fermions are neutral, we expect environ- 
mental coupling is weak. In addition, the lowest order 
local interaction term within the Majorana fermions is 
/ y(x)d x 'y(x)d^ r y(x)d^j(x) , which involves spatial deriva- 
tives at sixth order and will be strongly suppressed at low 
temperature. Thus there is reason to expect the low tem- 
perature dephasing rate for the Majorana fermion edge 
states will be smaller than that of ordinary electrons. 




FIG. 2: (a) A point contact for neutral Majorana fermions 
characterized by reflection and transmission amplitudes t and 
r formed by a junction between two superconductors. Each 
superconductor is connected to a source and drain by chiral 
fermions at a magnetic domain wall, allowing t to be measured 
with charge transport, (b) Zero energy transmission of the 
point contact as a function of phase <j> for different coupling 
strengths. The insets indicate the limits of a weakly coupled 
point contact (right) and a long line junction (left). 



We next study the transmission of Majorana fermions 
across a Josephson junction between two superconduc- 
tors, shown in Fig. 2a. The junction plays the role of 
a point contact for Majorana fermions and can be char- 
acterized by a scattering matrix relating incoming and 
outgoing Majorana modes, 7° ut = Sf^E)^™. Each su- 
perconductor is connected to a source and drain by chiral 
electron modes at magnetic domain walls. An incident 
electron from Si splits into two Majorana modes. One of 
the two is scattered by the junction, and has a probably 
amplitude t = Sfi of being transmitted and recombining 
with its partner before going to D± . Following the previ- 
ous procedure, we calculate the scattering matrix relating 
an incident fermion at Si to an outgoing fermion at D\ 
to obtain the current flowing to D\ when Si is at voltage 
V and the other leads are grounded. 
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where 6(E) is the same as in (j4j. At E = particle- hole 
symmetry constrains S pc to be a real 0(2) matrix de- 
scribing the transmission t = cos S and reflection r = sin 5 
such that 7$ mt + i^ ut = e l5 (^ + i'yf). The zero bias, 
zero temperature conductance, G = Idx/Vsx = te 2 /h di- 
rectly measures the transmission of the neutral Majorana 
fermions at the junction. 

The transmission amplitude t can be controlled by ad- 
justing the phase difference <j> of the Josephson junction. 
t(4>) depends on the geometry of the junction. We con- 
sider a simple model, 

H = ( 11 , l2 )[~iv M T z d x + X(x)cos(<t>/2)T y \{ lul2 ) T . (7) 

When X(x) — XS(x) and X/vm < 1, fl describes super- 
conductors weakly coupled by single electron tunneling 
at a point 0,[2l [2J|. When X(x) = A for x € [0, L] and 
otherwise, H becomes the low energy theory of a line 
junction^]. The transmission amplitude at E — in this 
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model is 

t(<t>) = l/cosh[Ccos(^/2)], (8) 

with £ = J <1xX(x)/2vm- Fig. 2b shows t(4>) for different 
values of C- At (j> — it, the transmission is perfect. This is 
guaranteed by gauge invariance. When <fi — > <f> + 2ir one 
of the Majorana edge modes changes signQ so r((j>) = 
—r((j) + 2tt). Thus, r(cp) — and t(<fi) = 1 for some <fi € 
[0, 2ir]. For a symmetric junction this occurs at (f> = 7r. 

For a weakly coupled point contact (Fig. 2b, right 
inset), £((/>) is energy- independent, but is only weakly de- 
pendent on <fi. For a long line junction, (Fig. 2b, left 
inset) t(4>) varies over a wide range of values between 
and 1, but has a very narrow peak 5<p ~ Kvj\{/ AqL. In 
addition, near the peak the transmission will be strongly 
energy dependent due to the small gap when <j) ~ tt. 
It is desirable to engineer the size and geometry of the 
Josephson junction in between these two limits, so that 
t(c/>) has a well defined peak which can be probed by the 
low temperature conductance. 

It is worthwhile to compare the superconducting point 
contact for Majorana fermions studied here with a point 
contact in the v = 5/2 quantum Hall effect. Our point 
contact is precisely equivalent to the neutral sector of 
the v = 5/2 point contact, which has been described in 
terms of the Ising boundary conformal field theory(25l|. 
For v = 5/2, however, the physics is dominated by the 
backscattering of charge e/4 quasiparticles, which is anal- 
ogous to quantum tunneling vortices across the supercon- 
ductor in our system. Since the superconducting phase 
is essentially a classical variable, this process is strongly 
suppressed in a superconducting point contact. Thus, 
unlike the v = 5/2 problem, vortex backscattering does 
not lead to a crossover to the weak tunneling limit. 

The recently discovered topological insulator 
Bi 2 Se3[2(| |26[, which has a large bulk gap ~ .35eV 
is a promising material to probe these states. Unlike 
Bii^Sbz, its surface states have a small Fermi surface 
that encloses a single Dirac point. Photoemission 
experiments reveal a Fermi velocity Kvp ~ .3eV nm 
and a Fermi energy /i ~ .3eV relative to the Dirac 
point. The current materials are unintentionally doped, 
with the bulk Fermi energy in the conduction band. If 
the material can be compensated cither by doping or 
gating, it is likely that the surface Fermi energy can be 
made much closer to the Dirac point. This is important 
because achieving the magnetic gapped state requires 
a field M > fi. Moreover, the k ■ p theory predicts 
that the Majorana velocity Vm is suppressed when 
Ao <C |U, reducing the temperature scale T$l required to 
observe the signature of Majorana fermions. Our model 
calculation gives Vm ~ uf(Ao/^) 2 . Assuming a super- 
conductor can be found that gives a proximity induced 
gap Ao ~ .lmeV, we require size L > Hvf/^-o ~ 3/zm. 



If /Li ~ lmeV and SL ~ l^m then Tsl ~ 30mK. T$l can 
be larger if the path difference 6L can be finely tuned. 

To conclude, we have proposed experiments to probe 
the interference and transmission of neutral Majorana 
fermions with charge transport. We hope they offer a 
first step towards the more ambitious goal@ of detecting 
the non-Abelian statistics of individual Majorana bound 
states and using them for quantum computation. 

In a recent preprint, Akhmerov, et al.(27j indepen- 
dently studied an interferometer similar to Fig. 1. We 
thank Carlo Beenakker for an insightful discussion. This 
work was supported by NSF grant DMR-0605066 and 
ACS PRF grant 44776-AC10. 
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